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ABSTRACT 

We investigate the spatial distribution of the baryonic and non-baryonic mass components 
in a sample of 66 virialized systems. We have used X-ray measurements to determine the de- 
projected temperature and density structure of the intergalactic medium and have employed 
these to map the underlying gravitational potential. In addition, we have measured the depro- 
jected spatial distribution of galaxy luminosity for a subset of this sample, spanning over 2 
decades in mass. With this combined X-ray/optical study we examine the scaling properties 
of the baryons and address the issue of mass-to-light (M/L) ratio in groups and clusters of 
galaxies. 

We measure a median mass-to-light ratio of 224 /170 (M/L)q in the rest frame Bj band, 
in good agreement with other measurements based on X-ray determined masses. There is no 
trend in M/L with X-ray temperature and no significant trend for mass to increase faster 
than luminosity: M oc Lg ( j 8±0 - 12 . This implied lack of significant variation in star formation 
efficiency suggests that gas cooling cannot be greatly enhanced in groups, unless it drops 
out to form baryonic dark matter. Correspondingly, our results indicate that non-gravitational 
heating must have played a significant role in establishing the observed departure from self- 
similarity in low mass systems. The median baryon fraction for our sample is 0.162 h 70 1 , 
which allows us to place an upper limit on the cosmological matter density, Q m ^ 0.27 hj^, 
in good agreement with the latest results from WMAP. 

We find evidence of a systematic trend towards higher central density concentration in 
the coolest haloes, indicative of an early formation epoch and consistent with hierarchical 
formation models. 
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1 INTRODUCTION 

As the material most amenable to detailed study, the properties of 
optically luminous matter have long been used to infer the distri- 
bution of mass in the Universe. Clusters of galaxies provide an ex- 
cellent laboratory for this purpose, since they are sufficiently large 
that their global properties reflect those of the Universe as a whole. 

Previous work has shown that both optical light and gas, 
in particular, are more spatiall y extended than the d ark matter 
within virialized systems iDavid. Jones & Fo rman 1995). Evidence 
of a monotonic rise in gas fraction with radius has also been 
found across a wide ra nge of mass scales in a very la rge sam - 
ple of virialized systems 1 Sander son et all2003l hereafterlPaper j). 
and al so in the recent analysis of Tcastillo-Morales & Sch indler 
(2003), which is in broad agreement with numerical simulations 
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lEke. Navarro & FrerdJ Il99ll iFrenk et alJll999h . This behaviour 
has been further confirmed in a high quality XMM-Newton obser- 
vation of a relaxed cluster (Pratt & Arnaud 2002), where the X- 
ray halo has been traced out almost to the virial radius, although 
Chandra observations of extremely relaxed lensing clusters have 
revealed a flat gas fraction profile in a small number of cases 
lAllen. Schmidt & Fabiani r2002 ) . However, the latter observations 
only probe the innermost ~l/3 of the halo: our data, though of 
poorer quality, extend beyond this region in many cases, and we 
benefit from averaging over a large ensemble of virialized systems. 

While the hot gas has been shown to depart systemati- 
cally from self-similarity under the influence of non-gravitational 
physics, the same cannot be said of the dark matter, which 
is not directly affected by such processes. N-body simula- 
tions indicate that the dark matter of virialized haloes should 
follow a universal profi le, across a wide range of scale 
sizes iNavarro. Frenk & White! [l995), apart from a mild trend 



© 2003 RAS 



2 A. J. R. Sanderson and T. J. Ponman 



in its concentration with mass jNavarro. Frenk & White 


1997 


Salvador-Sole, Solanes & Manrique 1998; Avila-Reese et al. 


1999 


Jina 2000). 



Although the contribution to the baryonic mass component 
from the individual galaxies in virialized systems is significantly 
less than that from the X-ray emitting intergalactic medium, the 
galaxy component is of great importance, none the less. As the 
ultimate endpoint of gas cooling, as well as the source of non- 
gravitational energy injection into the intergalactic medium (IGM), 
the behaviour of the stars is closely linked to the properties of the 
gas. In addition, stellar sources are responsible for the synthesis 
of the metals that contaminate the IGM, which are relea sed vi a 
supernova-driven winds (e.g. lFinoguenov. Arnaud & Davidl200lt) . 

However, the scaling properties of the stellar content are rather 
less well established than those of the IGM. In particular, there 
is some debate on the mass-to-light ratio in groups and clus- 
ters and how this compares with the global value. Recent stud- 
ies have reached differing conclusio ns about how this quantity 
scales with halo mass. For example. iGirardi et alJ 1200*21) report a 
trend towards increa s ed M l L in clusters, a result also favoured by 
Marin oni & Hudson! |2002). By contrast, measurements of M/L 
based on X-ray mass estimates have concluded that it is ro ughly 
universal in groups a nd clusters ICirimele. Nesci & Trevesel 1997k 
iHradeckv et all200d) . 

The approach of a combined X-ray and optical analysis of 
clusters and groups of galaxies is a powerful tool for the under- 
standing of mass distribution in these systems. The X-ray data yield 
information about the behaviour of the diffuse, hot intergalactic 
medium (IGM) and can be used to determine the underlying grav- 
itational potential structure, under the assumption of hydrostatic 
equilibrium. By contrast, optical measurements can be used to map 
the spatial properties of the baryons resident in luminous matter, 
and he nce dedu ce their contribution to the total mass budget. 

In lPaper we presented a detailed X-ray analysis of the IGM 
in a large sample of 66 virialized systems, allowing us to recon- 
struct the properties of the IGM. We have assembled our sample 
from three existing X-ray studies of groups and clusters, to which 
we have added a small number of cool groups. To each system, we 
have fitted analytical profiles to parametrize both the gas density 
and temperature as a function of radius. This allows us to place all 
the X-ray data on a unified footing, giving us the freedom to ex- 
trapolate the gas properties to arbitrary radius. As with detailed 
X-ray analyses, previous combined X-ray/optical studies of this 
nature have been restricted to relatively s mall sample sizes (e.g. 
Cirim ele"et all 19971 : IHradeckv et all2000l) . In this paper, we build 
on our modelling of the IGM, by incorporating the spatial distri- 
bution of stellar material, which also allows us to determine the 
distribution of dark matter. Despite only covering half our original 
sample (some 32 groups and clusters), our optical sample is well 
suited to the study of the scaling properties of these systems, since 
we retain good coverage across a wide range of system masses. 

Throughout this paper we adopt the following cosmological 
parameters; Ho = 70 km s _1 Mpc -1 and go = 0. All quoted er- 
rors are la on one parameter, unless otherwise stated. 



2 3D GALAXY DENSITY CALCULATION 

The task of mapping the stellar mass in virialized systems is com- 
plicated by the discrete nature of the individual galaxies - unlike 
the evenly distributed gas, which smoothly traces the underlying 
gravitational potential. In addition, the contrast against the back- 



ground of a cluster or group of galaxies is much lower at optical 
wavelengths, scaling only linearly with number density, rather than 
in proportion to p gas 2 in the case of X-ray emission. These two 
difficulties are further compounded by the large angular size sub- 
tended by a typical group or cluster, which may exceed one degree 
on the sky - significantly larger than can be viewed by many cur- 
rent generation wide-field CCD cameras. Correspondingly, there 
have been comparatively few detailed photomet ric studies of r ich 
clusters made with CCDs (e.g. Carlbe rg. Yee & E llingson 1997). 

One way to avoid some of the above problems is to quanti- 
tatively establish cluster membership with redshift measurements 
of individual galaxies in its vicinity. Although this avoids the need 
to estimate the background contribution, such an advantage comes 
only at the price of long observing times, in order to obtain high- 
quality spectra for a large fraction of the galaxy members. For this 
reas on such studies are restricted to a sma ll number of clusters 
(e.g.lKoranvi & GelleJ20oil2000l;lMohr et alll99rilGirardi etail 
ll995tlFabricant. Kent & K urtzll 19891) . Given these limitations, we 
have taken a different approach to the issue of measuring the spa- 
tial distribution of galaxies and then estimating the total luminosity 
of the whole system. We outline below our method, which is based 
on widely-available digitized photographic plate data derived from 
all-sky surveys. 

The process can be separated into two stages. Firstly, determi- 
nation of the spatial distribution of the galaxy number density (see 
section l2~3l and, secondly, the calculation of the normalization of 
this density profile (see section l3~2t . Although the normalization 
can, of course, be calculated in the spatial fitting, this makes no al- 
lowance for the contribution to the total luminosity from galaxies 
too faint to be observed. To estimate this contribution it is neces- 
sary to determine the luminosity function of the observed galaxies 
and extrapolate this down to some limiting magnitude, in order to 
correct for the missing light. Whilst it is possible to do this using 
photographic plate based data, the measured magnitudes require 
careful calibration, due to subtle variations in the sensitivity of the 
photographic emulsion both within and between plates. Therefore, 
we have taken integrated luminosity values from a small number 
of sources in the literature and used these to infer a 3-dimensional 
normalization, in combination with our fitted profile parameters. 

We assume a distribution of optical light described by an NFW 
profile jNavarro et alJl9 9S). i.e. 



where x — r/r s and r s is a characteristic scale radius. This func- 
tion rises from a logarithmic slope of -3 at large radii to a central 
cusp (p oc r -1 ) with the transition between the two regimes oc- 
curring around r s . Such a parametrization is advantageous as it has 
only one free parameter (since the normalization is constrained by 
our maximum likelihood fitting method - see section l2~3l . which 
helps stabilize the fitting, given the relatively large background 
level in the case of the plate data. 

We make the assumption of spherical symmetry in the depro- 
jection of the optical light distribution, as we did for the X-ray anal- 
ysis in lPaper j Since our sample has been selected on the basis of a 
relaxed X-ray morphology, it is reasonable to expect a correspond- 
ing degree of regularity in the galaxy distribution. This also im- 
proves the quality of the fit, which might otherwise be degraded by 
the presence of significant substructure or bi-modality. However, 
despite this, it was not possible to obtain a satisfactory surface den- 
sity fit to the photographic plate data for a small number of cool 
groups; their treatment is described in section l2~2l below. 
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Name 


RA 


Dec. 


z 


'■pa 


R200 


< 


i 

iJ 


M/L 


Reference c 




(J2000) 


(J2000) 




(keV) 


(arcmin) 


(arcmin) 


do 11 Ln • 0) 


/ AV 1S,J 




NGC 4325 


185.825 


10.622 


0.0252 


0.90 


22.4 


2.6 


1 99 + 0.6 

>■•• "-a. 6 


293+160 
zyo -160 


G02 (APS) 


NGC 5846 


226.385 


1.696 


0.0058 


1.18 


95.5 


39.3 


1 «n+ ' 95 

1 - 80 -0.9S 


9 04 +107 
zu *-107 


H02 


HCG 62 


193.284 


-9.224 


0.0137 


1.48 


33.5 


12.7* 


2.5911;! 


78+ 40 
'°-40 


G02 (APS) 


NGC 5129 


201.150 


13.928 


0.0233 


1.54 


20.5 


7.8 


2.5611J 


09+48 


G02 (APS) 


NGC 2563 


125.102 


21.096 


0.0163 


1.61 


31.4 


12.1 


1.76l°;» 


161+g 


H02 


Abell 262 


28.191 


36.157 


0.0163 


2.03 


50.4 


47.0* 


9.7llJ;» 


1181" 


G02 (APS) 


Abell 194 


21.460 


-1.365 


0.0180 


2.07 


52.2 


20.4 


5.16+1° 


090+2OO 
ozu -200 


G02 (COSMOS) 


MKW4 


180.990 


1.888 


0.0200 


2.08 


35.5 


40.3 


1 12+ ' 6 


glo+310 
DiJ -310 


G02 (APS) 


MKW4S 


181.647 


28.180 


0.0283 


2.46 


28.9 


11.9 


3.37 + }'2 


O9Q+130 


H00 


NGC 6338 


258.825 


57.400 


0.0282 


2.64 


26.4 


5.9 


14.01™ 


59 +60 
Oy -50 


G02 (APS) 


Abell 539 


79.134 


6.442 


0.0288 


2.87 


38.0 


5.9 


6.i7ll:l 




GOO 


AWM4 


241.238 


23.946 


0.0318 


2.96 


40.2 


11.5 


4.75+tl 


8861^ 


G02 (APS) 


Abell 1060 


159.169 


-27.521 


0.0124 


3.31 


104.8 


15.0* 


18.9HJ 


244 +i20 
z **-120 


G02 (COSMOS) 


Abell 2634 


354.615 


27.022 


0.0309 


3.45 


32.1 


12.0 


25.411*1 


76 +40 


G02 (APS) 


Abell 2052 


229.176 


7.002 


0.0353 


3.45 


36.6 


2.3 


6.7o!l;£ 


5 8 q+240 


H00 


Abell 2199 


247.165 


39.550 


0.0299 


3.93 


34.3 


31.2* 


9.98124 


211+ 60 


H00 


Abell 2063 


230.757 


8.580 


0.0355 


4.00 


35.7 


2.3 


7.0811;; 


541+WO 
0211 -140 


H00 


AWM7 


43.634 


41.586 


0.0172 


4.02 


105.2 


4.2* 


7-68+1* 


1614 +970 


GOO 


Abell 3391 


96.608 


-53.678 


0.0536 


5.39 


26.9 


4.3 


oc 1 +12.6 
zo.±_ 12 g 


9-.O + 120 
zlo -120 


G02 (COSMOS) 


Abell 2670 


358.564 


-10.408 


0.0759 


5.64 


19.6 


5.9 


20 7+ 10 - 3 


240+130 
z * a -130 


G02 (COSMOS) 


Abell 119 


14.054 


-1.235 


0.0444 


6.08 


33.5 


5.5 


2fi Q+ 13 - 6 
zo.y_ 13 4 


9 1S +120 
zi8 -120 


G02 (COSMOS) 


Abell 496 


68.397 


-13.246 


0.0331 


6.11 


39.2 


20.0* 


1Q 1+6.6 


091+I7O 
,5Z1 -170 


G02 (COSMOS) 


Abell 3558 


201.991 


-31.488 


0.0477 


6.28 


28.7 


18.2 


A7.7±f 3 % 


qq+50 


G02 (COSMOS) 


Abell 3571 


206.867 


-32.854 


0.0397 


7.31 


39.8 


17.9 


46.8li;t 


igi+80 


G02 (COSMOS) 


Abell 2218 


248.970 


66.214 


0.1710 


8.28 


11.3 


10.4* 


47 5+ 5 ' 9 


167+ 38 


S96 


Abell 1795 


207.218 


26.598 


0.0622 


8.54 


28.3 


6.4* 


9 7 Q+13.9 


oon+240 
■3^U_ 240 


G02 (COSMOS) 


Abell 2256 


256.010 


78.632 


0.0581 


8.62 


27.3 


7.3 


09 9 + 19.6 


1 75+90 
1 '°-90 


GOO 


Abell 85 


10.453 


-9.318 


0.0521 


8.64 


28.1 


11.6 


no n+12.1 
zo.»_ 12 Q 
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zzy -120 


G02 (COSMOS) 


Abell 3266 


67.856 


-61.417 


0.0545 


9.53 


29.7 


12.0 


29.3ll 4 4 ? 


26 i+140 
ZD1 -140 


G02 (COSMOS) 


Abell 2029 


227.729 


5.720 


0.0766 


9.80 


26.4 


4.6 


86.3+^? 


1 55+8O 


G02 (APS) 


Abell 478 


63.359 


10.466 


0.0882 


10.95 


17.8 


12.8* 


28.1 +7 7 l 


210+ 160 
ZiU -160 


H00 


Abell 2142 


239.592 


27.233 


0.0894 


11.16 


22.7 


8.7 


5o.oif 5 :° 


25Q+170 
Z3U -170 


G02 (APS) 



Table 1. Some basic pro p erties of the 32 obje cts in the optical sample, listed in order of increasing temperature. Positions and redshifts are taken from 
Ebelin getalJll996l l 1998); Ponman et al. 1 1996) and NED. Columns 7-9 are data as determined in this work, except for those values of r s marked with a *, 
which were taken from Lloyd-Davies i200l|). Note: values are for Hq = 70 km s —1 Mpc^ 1 . All errors are 68% confidence. "The cooling-flow corrected, 



emission-weighted temperature of the system within 0.3i?2O0, as determined in Paper 1. The NFW scale radius of the galaxy distribution (see equation 



_>0Q , : 

Refe rence for the lumino sity data, G02 = lGirardi et afl<2002l) : GOO = lGirardi et alj|200fj) : H00 = lHradeckv et alj<200fj) : H02 = lHelsdon & Ponmanl 
S96 = ISauires et ail Il996ll . 
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2003); 



2.1 APM data 

The Automatic Plate Measuring (APM) machine source catalogue 1 
is a digitized catalogue based on the 1st generation Palomar Obser- 
vatory Sky Survey (POSS I) and United Kingdom Schmidt Tele- 
scope (UKST) photographic plates. The field of view covered is 
extremely large - 6.2 and 5.8 degrees on a side for the plates from 
POSS I and UKST, respectively - and the whole sky has been ob- 
served in this way. As such, it is ideal for studying large scale struc- 
ture at optical wavelengths, although the magnitude information 
it provides can be rather unreliable, owing to saturation and other 
problems with photographic emulsions. However, we only exploit 
the positional data available - which are very reliable - since we 
determine total luminosities for our groups and clusters from liter- 
ature measurements. 

To extract source positions for our fitting, we have applied 
very conservative selection criteria to allow for problems with 



http://www.ast.cam.ac.uk/~apmcat/ 



source identification. Although detection algorithms have been 
used to distinguish between stellar and non-stellar sources, this 
classification is often unreliable and can cause stars to be identified 
as galaxies (e.g. Caretta, Maia & Willmer 2000) and vice versa: in a 
recent analysis of data from the Minnesota Automated Plate Scan- 
ner (APS) digitized survey of POSS I plates. iRines et al] 1 2002) 
found that ~3.3 per cent of objects identified as stars were actu- 
ally galaxies, in the vicinity of the cluster Abell 2199. This con- 
fusion is exacerbated for faint galaxies, and consequently we only 
neglect sources classified as stars if they are brighter than = 15, 
otherwise all detected sources are used. This avoids any potential 
systematic bias in the fitting between clusters, at the expense of a 
significantly increased background level. All objects classified as 
noise features are excluded. 

2.2 NED data 

The problems associated with optical analysis are exacerbated at 
the scale of groups of galaxies, which present even lower contrast 
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against the background, owing to their smaller volume. However, 
many well-studied groups - including most of those in our sample 
- are sufficiently close so as to have measured recession veloci- 
ties available in the literature. In such cases it is therefore possible 
to establish unambiguous membership of the group, thus eliminat- 
ing the contamination from fore- and background objects. How- 
ever, the numbers of such confirmed members are generally rather 
small and tend to be predominantly limited to the inner regions 
of the halo. The impact this has is mitigated to some extent by 
our chosen parametrization of galaxy density (equation 0, which 
has only a single free parameter (the scale radius, r s ), since sparse 
data obviously provide weaker constraints on spatial fitting (e.g. 
lGirardietalll995h . 

It was not possible to derive NFW scale radii from the APM 
data for four of our groups (NGC 2563, NGC 5129, NGC 5846 
& NGC 6338) and so we took galaxy positions from the NASA 
Extragalactic Database (NED) for these systems. We used only 
those galaxies which were identified as being group members in 
the database. It was also necessary to do this for Abell 539, since 
no POSS I or UKST data were available for this poor cluster. There 
is potential for bias when using NED data, since the completeness 
properties of the database are unknown - although we are only 
concerned with the spatial completeness, since the luminosities for 
these systems are taken from other sources. In particular, it is pos- 
sible that there may be a dearth of galaxies catalogued in the outer 
regions of these groups, which may introduce a spurious central 
surface density enhancement. However, we note that any such bias 
is reduced by the fact that we have used the radius of the outermost 
galaxy as the cut-off radius in the fitting, rather than the X-ray de- 
termined virial radius used to select the galaxies from NED. 

2.3 Surface Density Fitting 

As we are interested in the 3-dimensional properties of the galaxy 
optical light, we have chosen to parametrize the space density dis- 
tribution in three dimensions and then numerically project this, 
for comparison with the surface density of galaxies. The same ap - 
proach was used in the deprojection of the X-ray data (see Paper I). 
This method is advantageous, since it directly yields the required 
parameters, without the need for inversion. For the analysis of the 
APM data, a constant background term was included as a free pa- 
rameter, to allow for a variation in the number of fore- or back- 
ground sources caused, for example, by a differences in large-scale 
structure along the line-of-sight. It is important to fit this separately 
for each cluster, since our background (incorporating a stellar com- 
ponent) can vary significantly with position on the sky. The centroid 
of the model was fixed at the cluster position as listed in NED, in or- 
der to stabilize the fitting. Although the optical centroid can be dis- 
placed with respect to the X-ray centroid, we note that lGirardi et alj 
1 2000) find a typical error in Lb within R v of just 5 per cent, when 
recomputing their integrated optical luminosities using an X-ray 
rather than optical centroid. 

The 3-dimensional density profile was evaluated in a series 
of spherical shells, by numerically integrating equation Q between 
their radial boundaries. For each spherical shell there is a corre- 
sponding annulus, with identical radial bounds. The surface density 
in each annulus is calculated by summing the contributions made to 
it by its corresponding shell and all those lying outside it. Hence the 
contribution to the surface density from projection along the line of 
sight is fully accounted for. 

We hav e chosen t o perform an u nweighted fit to the data (as 
also used bv lHradeckv et all2000l and lKoranvi & Gellej20oj for 



example), so that each galaxy is treated as a single point at a given 
position. This approach is advantageous since it avoids the need to 
use plate-measured magnitudes, which require careful calibration 
as mentioned above. By treating each galaxy identically we are im- 
plicitly assuming that there is no luminosity segregation, i.e. that 
the luminosity function is everywhere the same within the cluster. 
This is a reasonable assumption, since any segregation is generally 
limited to only the very brightest gal axies - which tend to be lo - 
cated more centrally in clusters ( Adami, Katgert & Biviano 1998). 

The use of an unweighted fit has important consequences for 
those clusters that possess a central cD galaxy. On the one hand, 
our method limits the potential for an extremely bright, centrally 
located galaxy to skew the fit, by overemphasising the central lu- 
minosity density. On the other hand, in so doing, there is a tendency 
to slightly overestimate the luminosity in the outer regions. This is 
due to our method for determining the luminosity normalization, 
which does include a contribution from any central cD object. As 
a result, the extra luminosity of such a galaxy is, in effect, 'spread 
out' over the whole profile, leading to an artificial excess in the 
outer regions. However, the effect this has on our total luminosity 
estimates is small, given the generally large sizes of the apertures 
used to measure the luminosity values to which we normalize our 
density profiles (see section l3~2l . In any case, the contribution to 
the luminosity function from cD galaxies is fully accounted for by 
this same method for determining the total cluster luminosity. 

Although some previous studies of galaxy distributions in 
cluster s have fitted analytical profiles wi th cores to the data (e.g. 
iGirardi et al1ll995t Icirimele et alJll997l) . we favour the cusped 
NFW profile (equation^. Not only does it have one fewer fitted 
parameter - which helps stabilize the fitting, as previously men- 
tioned - it also pr ovides a more accur ate description of the dark 
matter distribution iNavarro et alll995h . Simulations indicate that 
the galaxies trace the dark matter more closely than the gas does, 
with only a modest bias iDave. Katz & Weinber J2002I) : they may 
even be more centrally concentrated, owing to the effects of dynam- 
ical fr iction iM etzIer & Evrardfi997l) . We note that lAdami et alJ 
1 1998) have shown that only the surface density of faint galaxies 
shows a significant preference for a distribution with a core rather 
than a cusp. Since our analysis is based on photographic plate ob- 
servations, rather than the higher quality CCD data used in their 
work, we can therefore expect to be relatively insensitive to this 
effect. 

Previous studies offer support for the presen ce of a cusp in the 
lumino sity distribution. For example, Oege rle. Jewison & Hoessell 
1 1987) measure a power law slope of a pproximately -1 in the cen- 
tres of 3 Abell clusters; Merritt & Trembla^ J 1994) find evidence 
of a central power law dependence of galaxy space density in 
the core of the Coma cluster, using a non-parametric algorithm. 
iBeers & TonrvJ ll986h have demonstrated that imperfect centring 
can 'erase' the presence of a real cusp and artificially introduce a 
constant density core. Using two methods of determining the clus- 
ter centre, they find that the majority of the 48 clusters in their sam- 
ple posses a central cusp in projected galaxy density. The inner re- 
gions of these clusters are well characterized by a power law slope 
of-1. 

We use th e maximum likelihood fitting procedure described in 
Sarazi 3 <198Cl) . which we now briefly outline. There are three pa- 
rameters to be determined - the scale radius, r s , the surface density 
normalization, Eo, and background surface density, Eb g (since the 
centroid is fixed - see above). However, only two of these parame- 
ters are independent, therefore we require an additional constraint, 
otherwise the maximum likelihood solution is unbounded. We ap- 
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ply the condition that the model reproduce the observed number of 
galaxies within the fitting region. In the standard method of lSarazinl 
1 198C), r s and So are left as free parameters and Eb g is calculated 
so that the 'area' under the model + background equals the total 
number of galaxies within the fit region. However, in the case of 
our NED data Eb g is known a priori (i.e. it is zero, since all the 
galaxies are confirmed group members], so we use the alternative 
method detailed in section IV of lSarazinl (l980). Here, Eb g and r s 
are left as fitted parameters, while Eq is instead calculated so as to 
reproduce the observed galaxy count. 

As a result of our numerical projection technique, Eo is not 
actually a direct model parameter. However it is trivially related to 
po, the NFW central density normalization from equation Q We 
determine the appropriate value of po by setting it initially to unity 
and calculating the 'area' under the projected data, excluding the 
parts of any annuli lying outside the fit region, i.e. summing the 
product of the surface density in each annulus and its area. The 
ratio of the observed number of galaxies in the fit region to this 
number then yields po- 

Fitting was performed using the gradient-based MIGRAD 
metho d in the MINUIT minimization library from CERN Jjamesl 
1 19981) and errors on the parameters were found with MINOS, from 
the same package. Errors were determined from the increase in 
the Cash statistic of one, since differences between values obtained 
from the same data set are % 2 -distributed. 



3 OPTICAL LUMINOSITY CALCULATION 

Since we have chosen to avoid using magnitude information in the 
calculation of the galaxy surface density distribution, we require 
an alternative means of deriving a normalization value for our 3- 
dimensional luminosity profiles. To do this we have taken inte- 
grated luminosity measurements in fixed metric apertures from the 
literature, which allow us to infer a central luminosity density. The 
details of this procedure are described in section l3~2l below. 

The majority o f these values are taken from the sample of 
which incorporates most of the sample of 



Girardietal. (2002) 
Girardi et al. 1 200G). This sample comprises data from the APS and 
COSMOS/UKST surveys - we have selected luminosity estimates 
from the latter for preference, where overlaps occurred. The only 
exception to this is HCG 62, for which the COSMOS aperture lu- 
minosity quoted in Girardi et al. 1 2002) is a factor of 2 larger than 
the APS luminosity; for this group we have chosen the APS mea- 
surement. We note that the APS luminosity leads to a mass-to-light 
ratio sec tion 14. 2\ which is more c onsistent with the simple esti- 
mate of iPonman & Bertrarnl ll993l) . Data for three other clusters 
(Abell 539, Abell 2256 and AWM 7) were taken from lGirardi et alJ 
(2000). These latter systems have two different luminosity esti- 
mates based on literature magnitudes - we chose the 'M(red)' sam- 
ple, since this provided the largest overall number of galaxies de- 
tected for all three clusters. 

A further five systems are covered by the study of 
iHradeckv et alJ l2000S) . based on data from the second generation 
Palomar Observatory Sky Survey (POSS II), calibrated with CCD 
photometry. Two grou ps (NGC 2 563 and NGC 5846) have lumi- 
nosity estimates from Helsdon & Ponman <2003l) within an aper- 
ture defined by the virial radius as calculated in Helsdon & Ponman 
2000). Finally, we used the luminosity quoted in ISauires et alJ 
1996) for the rich cluster Abell 2218. The references for each sys- 



3.1 Conversion between bands 

We have adopted the Bj photometric band as our standard reference 
frame , since this was used bv lGirardietaL[ j2002 ) and Girar dlet alJ 
(2000), and we take B Q j = 5.33 jGirardi et alJl2002l) . However" 
our other literature sources have used different bands and so we 
have converted these values into the Bj band, under the assumption 
that the majority of the light originates in early-type galaxies. We 
have used the following relations to perform the conversions. 
To convert from the B band to the Bj band, we assume 

B Q = 5.48 lGirardietalf2 000) 



Bj = B — 0.28f B — V) iBlair & Gilmorel 19821) 

B — V = 0.9, for early-type galaxies jGirardi et all200d) 

Bj = B - 0.252, 



which leads to 

Lb, Lt> 
- — — = 1.10—5- 



J B,0 



(2) 



To convert from the V band to the Bj band, we assume 
Vq = 4.82 1 Alien! 19731) 

=> Bj = V + 1.152, from above (for early-type galaxies) 
which leads to 

' - 0.554—^. 



l b q . ^V,0 

For comparison with R band luminosities, we assume 
R Q = 4.28 lAUenll973l) 



(3) 



V — R = 0.55 for early-type galaxies l Girar di et all200d) 
=> Bj = R + 1.702, from above, 

which leads to 

Lb, Lt> 
- — — = 0.547-—^-. 



"R,0 



(4) 



3.2 Determination of luminosity normalization 

With the exception of the IHradeckv et all {2000) data, all our lit- 
erature values of luminosity are simple estimates within a fixed 
aperture on the sky. As such, these measurements represent the in- 
tegrated luminosity within a cylinder, defined by the aperture mul- 
tiplied by the diameter of the cluster. It was therefore necessary 
to allow for the projected contribution to this value from luminos- 
ity at large radii, when inferring values for the 3-dimensional nor- 
malization. We accomplished this in the following way, using our 
numerical projection technique described previously. The normal- 
ization was set to unity and the NFW profile was integrated out 
to a large radius, to insure convergence (taken to be 500 arcmin), 
in a series of fixed-width spherical shells. The projected emission 
from these shells was then computed, in the corresponding series 
of annuli. The total luminosity within the aperture radius was then 
calculated, by simply summing the contributions from the annuli 
within that radius. Where the aperture radius lay within an annulus, 
the contribution from this annulus was determined by linear inter- 
polation, since the annulus width is small. The appropriate normal- 
ization factor was found by dividing the measured luminosity value 
by that calculated for a normalization of unity. 



tern are listed in TableQ together with some key properties. 



For the five clusters analysed by Hradec kv et alJ {2000), a 
slightly different approach was taken, since the authors corrected 
their aperture luminosities for the effects of projection. Working in 
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three dimensions, we set our normalization to unity and simply in- 
tegrated our density profile from r = to a radius of lh.50 Mpc 
(used for all their quoted luminosities). The ratio of the aperture lu- 
minosity to the integrated value then yields the central luminosity 
density that we require. 

The errors on our quoted luminosities are derived from the 
errors on the aperture values. Following Girardi et alJ 120021) . we 
assume a la error of 50 per cent on our total Lbj values based 
on their data. The same frac tional error was assumed for the val- 
ues from lHelsdon & Ponmanl <2003l). based on their own error esti- 
mates. The values for the Hradec kv et al. (2000) data have smaller 
errors, as does the lSauires et al]ll996l) luminosity for Abell 2218, 
since these are based on higher quality POSS II and CCD obser- 
vations, respectively, which are more accurate than the majority of 
our photographic plate data. 

Owing to our choice of an NFW profile for the fitting, coupled 
with the bias towards optically rather than X-ray selected clusters 
in previous studies, we are extremely restricted in the number of 
direct comparisons which can be made between our spatial fitting 
results and those available in the literature. However, we note that 
lAdami et alJll998Tl found a scale radius of 3.7 ± 0.6 arcmin for the 
cluster Abell 1 19, which is reasonably consistent with our value of 
5.5 arcmin (see Table0. This change in r s would give a ~4 per 
cent change in the total optical luminosity within R 200, based on 
the aperture luminosity for this cluster as measured by Girar di et alJ 
(2002). The impact of uncertainties in the spatial fitting parame- 
ters is limited by the fact that the apertures from which the optical 
luminosities we use for calibration have been taken, are close to 
-R200: the median difference between these aperture radii and our 
X-ray determined values of -R 200 is ~50 per cent . In the case of 
Abell 119, our correction to the Girardi et al. 12002) aperture lumi- 
n osity amounts to only ~10 per cent. For Abell 2256, the analysis 
of lOegerle et"al] l^Sl) provides support for our findings, measur- 
ing an inner logarithmic slope of —0.98 ± 0.02 for the galaxy den- 
sity distribution in this cluster. Moreover, their binned radial profile 
is consistent with an outer slope of -3 and a break radius similar to 
our value of 7.3 arcmin. 



4 RESULTS: OPTICAL PROPERTIES 

We convert optical light directly into stellar mass, assuming a mass- 
to-ligh t ratio for ea rly-type galaxies of (M/L)b = 7 hro Solar 
(Pizzella et al. 1997). A similar value of (8.3±0.35) h ro was found 
by Ivan der Marell ll99lh . Using Bq — Bq j = 0.15, from sec- 
tion o gives (M/L)bj = 6.1 hro{M / 'Lbj)0. To infer the den- 
sity profile of dark matter, we subtract the gas and stellar mass from 
the total mass (as determined from the X-ray data). Consequently, 
we only have information on the dark matter and stellar properties 
for the optical sub-sample of 32 groups and clusters. 

Relationships between integrated properties provide an impor- 
tant tool for investigating the similarity between haloes across a 
range of masses. Although quantities such as luminosity are di- 
rectly observable, some allowance must be made for the effects 
of projection if fair comparisons are to be made between differ- 
ent systems. Correspondingly, all our quoted luminosity values (see 
Table □ are derived from integrating the 3 -dimensional light pro- 
files, using the appropriate normalization calculated in section lX2l 
to provide a more sensitive probe of scaling properties. Where we 
have compared our results with those from the literature quoted in 
different photometric bands, we apply the correction factors given 
in section l3~T1 to determine equivalent values in the B s band. 




Total Mass (Solar) 



Figure 1. Total gravitating mass as a function of Lbj luminosity (both 
within -R200). f° r the optical sample. The line shows the best-fitting power- 
law, which has a logarithmic slope of 0.96 ± 0.10. 



4.1 Optical luminosity 

A useful test of the scaling properties of the stellar distribution 
is the relationship between optical luminosity and mass. Fig. Q 
shows the Lbj luminosity within -R200 for the optical sample plot- 
ted against the total mass within R200, as calculated in IPaper 5 
We performed an orthogonal distance regression to fit a straight 
line to the data in log space, using the ODRPACK software package 
I Bogg s et alll989l|l992l) . to take account of errors in both X and Y 
directions. The best-fitting relation is log Lbj = (— 1.36±1.46) + 
(0.96 ± 0.10) x log M and the scatter in the data is 1.21 times that 
expected from the statistical errors alone. The logarithmic slope of 
this relation is flatter than, but consistent with self-similarity (i.e. 
Lbj oc M), at the la level - this reflects the slight excess in stellar 
density in the coolest systems referred to in section l4~3l Perform- 
ing the regression with the axes reversed leads to a scaling between 
light and mass of M oc L^f ±0 - 12 . 

For comparison, we have also plotted Lbj against the 
emission-weighted temperature, as calculated in lPaper fl which we 
show in Fig.|2| Using the same regression technique, the best-fitting 
relation is log L B j = (11.02 ± 0.01) + (1.62 ± 0.14) log kT, and 
the scatter in the data is 0.96 times that expected from the statistical 
errors alone. Furthermore, there is no evidence of any steepening in 
the relation or any other systematic deviation from a simple power 
law, indicating that this provides an accurate description of the data. 
The logarithmic slope of this relation is steeper than, but consistent 
with self-similarity (i.e. Lbj oc T ' 5 , assuming that light traces 
mass), at the la level. However, this is due to the effect of the 
M — Tx relation slope - which is itself significantly steeper than 
the self-similar relation ( Paper I). 

The scatter about the Lbj — Tx relation is rather less than that 
about the Al — Tx relation, especially considering that our mass 
errors are quite conservative I Paper I). This is puzzling, but may 
reflect the fact that the optical luminosity is dominated by the stel- 
lar contribution from the inner regions, where the gas temperature 
is also more heavily weighted. This would produce a tighter corre- 
lation between Lbj and kT than between Lbj and M, since the 
total mass is less sensitive to the contribution from the inner regions 
of the halo. 
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Figure 2. Total Lbj luminosity within -R200 as a function of system tem- 
perature for the optical sample. The line shows the best-fitting power-law, 
which has a logarithmic slope of 1.62 ± 0.14. 
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Figure 3. Mass-to-light ratio within R200 (in the Bj band) as a function 
of system temperature for the optical sample. The solid line represents the 
best fitting power law, which has a logarithmic slope of —0.11 ± 0.18. 
The dotted line shows the median value for the whole sample, of 224 ft, 70 

(M/L B ,j) . 



4.2 Mass-to-light ratio 

Since light is easily observed, it is often used as a tracer of mass in 
the Universe, expressed as the ratio of total mass to optical lumi- 
nosity in some photometric band. Assuming that clusters of galax- 
ies are a fair representation of mass composition on large scales, 
the mass-to-light ratios in these systems can be used to estimate the 
total mass density of the Universe. It is therefore important to un- 
derstand the scaling properties of this quantity if such cosmological 
inferences are to be unbiased. 

Fig. [3] shows the variation in mass-to-light ratio with X-ray 
temperature, for the optical sample (see Table Q for the data). It 
can be seen that there is little evidence of any trend in the data; 
Kendall's K statistic indicates the significance of a negative corre- 
lation is 0.37a. The best-fitting power law has a logarithmic slope 
of —0.06 ± 0.17, which is fully consistent with no trend. The level 
of scatter about the best-fit is exactly that expected from the sta- 
tistical errors. We therefore conclude that a universal mass-to-light 



ratio in groups and clusters provides a good description of most of 
our data. 

It can be seen that one point in particular is very high 
on this relation: the poor cluster AWM 7 has (AZ/Lbj) ~ 
1600(M/Lb,j)q albeit with rather large errors. A large mass- 
to-ligh t ratio was also measured for this system by iKoranvi et alJ 
1 1998), who found a value of (650 ± 150) hi 00 (M/L)© for the 
R band, corresponding to (830 ± 190) h?o (M / Lb,j)q, consistent 
with our lcr lower bound of ~650. Exclusion of this point lowers 
the mean of our sample from 308±53 to 266±32 h^o (M/Lbj)©. 
We also note that the mass-to-light ratio for HCG 62 is rather low 
- only 78 ± 40, which is rather less than the value estimated in 
IPonman & Bertranll 19931) . As was mentioned in section[3] there is 
a factor of two d ifference between the aperture luminosities quoted 
in lGirardi et all §002), based on APS and COSMOS data for this 
system. Furthermore, we note that the optically-determined virial 
radius used for their aperture is twice as large as our X-ray mea- 
sured -R200, which means that we have to reduce the (Girardi et alJ 
( 2002) luminosity value, according to our fitted galaxy density pro- 
file. 

To minimize the bias caused by such discrepant measurements 
in our M/L data, we have evaluated a logarithmic mean value, of 
h 70 (M/Lbj)©. This is significantly lower than the ordi- 
nary, arithmetic mean value of 308 ± 53, but compares very well 
with the median value of 224 /170 (M/Lbj)©. This illustrates the 
effect that anomalously large M/L systems in particular can have, 
when determining a representative average. 

In a combined X-ray /optical study of 12 Abell clusters of 
galaxies, Cirim ele"et alJ 1 19971) measured a mean mass-to-light ra- 
tio equivalent to 346 /t 7o (M/Lb.i)^, which is consistent with our 
results. More recentlv. lFIradeckv et alJiToOOh have also determined 
M/L for a sample of 8 groups and clusters of galaxies, correcting 
their optical luminosities for projection effects, by assuming a spa- 
tial distribution of light identical to that they measured for the X-ray 
gas. Their mean value of M/L is equivalent to 310 ± 45 (M/L)q 
in the Bj band, which is in excellent agreement with our own arith- 
metic mean. Furthermore, their media n M/L of ^250 is in distin- 
guishable from our own median value. Ffrad eckvet all (2000) also 
conclude that M/L i s rough ly independent of mass, as was also 
found by iDavid e t al. ( 1995), albeit based on a rather small sam- 
ple in both cases. We note however, that their assumption that the 
optical light traces the gas mass is not supported by our data (see 
section l4~4l . although this is likely to have only a small effect on 
their results. 

The agreement with M/L as measured using optical mass es- 
timate s is generally l ess good, as pointed out by iHradeckv et all 
(2000). For example, Gir ardi et alJ {2000) quote a mean value of 
~ 175 /170 (M/Lbj)© for their sample of 105 clusters, from which 
many of our own luminosity estimates have been taken (since many 
of their results were incorporated in Gira rdi et alJ y002)). How- 
ever, the disagreement with our own mean value of 308 ± 53 hyo 
(M/Lb,j)© can be attributed to a difference in luminosity as well 
as mass: we have extended their quoted aperture Lbj values, ac- 
cording to our fitted galaxy density profiles, and evaluated an inte- 
grated luminosity within -R200 as derived from our X-ray data. On 
lower mass scales, opticall y-derived M / L measurements a re also 
lower than our own value. iRamella. Pisani & Gelled ^997) find a 
mean approximately equivalent to 151 ± 6O/170 {M/ Lbj )q (99% 
confidence) for a very large sample of groups in the Northern Cf A 
Redshift Survey, using virial mass estimates. 

The difference between X-ray and optically derived masses 
can also be seen in the scaling properties of M/L. While our re- 
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suits support a universal mass-to-light ratio in grou ps and clus- 
ters - consistent with the X-ray studies of iHradeckv et alJ fcOOOt) 
and lCirimele et alJ Jl997l^ - iGirardi etail <2002h find that mass in- 
creases more quickly than lu minosity, such that M oc Lg' 3 °' 03 . 
However, the optical study of Carlber g et alJ Il996l> indicates that 
M/L is universal in clusters, although their mean value, equiv- 
alent to 380 ± 70^70 (M/Lbj)q is somewhat higher than our 
own. Weak lensing prov ides the m ost robust determination of halo 
mass, and the analysis of Hoekstra et al. 1 2001 ) indicated that group 
mass-to-light ratios were lower than those of clusters, although 
their sample comprised only groups. Such a variati on of M/L 
with h alo mass is confirmed by the work of lBahcall & Comerford 
(2002), who incorporated lensing mass estimates in part of their 
sample. 

Semi-analytical models (SAMs) of galaxy formation gener- 
ally predict a s ignificant i ncrease in MIL with halo mass (e.g . 
iKauffman et alJll999tlBenson et alJl2000tlSomerville et alJl2001l). 
As an example of the typical behaviour observed, iBenson et alJ 
(2000) find that M/L reaches a minimum on mass scales of ~ 
1O 12 M0, increasing by a factor of 3 up to halo masses of 10 15 Mq . 
This corresponds to a logarithmic slope of 0.16, which is only just 
outside the la upper bound of 0.13 from our data in Fig. [3] al- 
though we only include groups and clusters in our optical sample. 
This tendency for M/L to increase away fro m a minimum on the 
scale o f a typical galaxy is also confirmed by M arinoni & Hudsonl 
(2002). Semi-analytical models also predic t that M/L in cluster s 
is significantly lower than the global value I Kauf fman et alll9"99l) . 
suggesting that even massive clusters are a biased estimator of the 
mean M/L in the Universe. That our findings do not provide a 
closer match to the M/L predictions of these models may reflect 
the over-cooling of gas in SAMs, which is a well-known problem 
that leads to the for mation of excessively bright cen tral galaxies, 
for example iKauffmann. White & Guiderdonil 19931) . 

4.3 Stellar Distribution 

The stellar mass density as a function of scaled radius, for the op- 
tical sample, is shown in Fig.|4| grouped into five temperature bins 
for clarity. There is a remarkably clear consistency in the shape and 
normalization of the profiles across the whole range of tempera- 
tures, which are approximately co-aligned. However, close inspec- 
tion of Fig. |4| reveals a slight excess in the stellar density of the 
coolest two temperature bands, compared to the more massive sys- 
tems. This points towards a possible weak trend in star formation 
efficiency with mass; we will revisit this issue in section l531 

4.4 Spatial variation of mass-to-light ratio 

The extent to which light traces mass within virialized haloes can 
be gauged by studying the mass-to-light ratio as a function of radius 
in these systems. Fig.|5]shows these profiles for the optical sample, 
grouped by temperature as before. It can be seen that there is evi- 
dence of decrease in this quantity with increasing radius, outside of 
~0.1T?200, in four of the five bands. Interestingly, the central range 
(2.9-A.6 keV) shows a clear increase in M/L with radius in this 
region, although this may be partly driven by the anomalous be- 
haviour of AWM 7, which exhibits a strongly rising mass-to-light 
profile, as also found by Koranvi et al. 1 1998). 

Outside of ~0. 1-R200 our results demonstrate that the distri- 
bution of optical light is more extended than that of the gravitating 
mass, as previously reporte d bvlDavid et alJll995l). I t is worth not- 
ing that the simulations of lMetzle^^vrara " ^?97l) indicate that 
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Figure 4. The variation of stellar matter density with scaled radius, for 
the optical sample, grouped by system temperature. The solid line repre- 
sents the coolest systems (0.3-2.0 keV), increasing in temperature through 
dashed (2.0-2.9 keV), dotted (2.9-4.6 keV), dot-dashed (4.6-8.5 keV) and 
finally dot-dot-dot-dashed (8.5-17 keV). 

the stellar distribution can be less extended than that of the dark 
matter, due to the effects of dynamical friction transferring energy 
away from the galaxies. 

We have chosen to truncate the radial scaling in Fig. [5] at 
0.037?200, rather than O.Olifeoo. as our optical luminosities are 
likely to be unreliable in the core of the halo. This is because our 
unweighted galaxy density fitting is insensitive to the excess lu- 
minosity associated with a central cD galaxy. However, we would 
therefore expect to underestimate the luminosity in the core, which 
would overestimate the mass-to-light ratio and hence flatten any 
drop in M/L in the centre. This suggests that the observed trend 
towards a decrease in the mass-to-light ratio < O.I-R200 is probably 
real. We also note that a similar decline in M/L towards the centre 
is see n in the central regions of most of the 12 clusters in the sam- 
ple offcirimele et al. 1 1997). A recent detailed Chandra analysis of 
Abell 2029 also reveals a sh arp decrease in M/L within 0. 1-R200 
Lewis. Buote & Stockel2003l) . 



5 RESULTS: DARK MATTER, GAS AND TOTAL MASS 
PROPERTIES 

Building on our analysis of the optical properties of our sub-sample 
of 32 groups and clusters, we now present equivalent information 
for the other mass components. Since we infer the dark matter dis- 
tribution using our knowledge of the stellar mass profile, we are 
limited to our 'optical sub-sample' in the analysis of its properties. 
However, the gas and total gravitating mass are determined from 
the X-ray data alone, and so here we are able to improve our statis- 
tics by using data from the full sample, rather than being limited to 
those systems for which we have optical data. 

5.1 Dark matter distribution 

The spatial variation of the density of dark matter can be seen in 
Fig. |6| for the optical sample. The similarity between the hottest 
four temperature bands is quite close, with only the coolest systems 
showing a deviation from the general trend, and then only within 
~0.17?2O0- The dark matter distribution is largely self-similar, as a 
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Figure 5. The variation of mass-to-light ratio (in the Bj photometric band) 
with scaled radius for the optical sample, grouped by system tempera- 
ture. The solid line represents the coolest systems (0.3-2.0 keV), increasing 
in temperature through dashed (2.0-2.9 keV), dotted (2.9-A.d keV), dot- 
dashed (4.6-8.5 keV) and finally dot-dot-dot-dashed (8.5-17 keV). 
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Figure 6. The variation of dark matter density with scaled radius, for 
the optical sample, grouped by system temperature. The solid line repre- 
sents the coolest systems (0.3-2.0 keV), increasing in temperature through 
dashed (2.0-2.9 keV), dotted (2.9-4.6 keV), dot-dashed (4.6-8.5 keV) and 
finally dot-dot-dot-dashed (8.5-17 keV). 

consequence of its insensitivity to the types of heating and/or cool- 
ing processes which can influence baryonic material. Despite this, 
it is clear that there is evidence of an enhanced central density in 
the core of the average density profiles of the coolest groups. We 
return to this unusual behaviour in more detail in section l5"4l 

5.2 Total density 

To provide some comparison with the density profiles of the sep- 
arate mass components we have plotted the profiles of integrated 
overdensity in Fig. Q i.e. the mean density within a given radius, 
normalized to the critical density of the Universe. Since the grav- 
itating mass profile is determined from the X-ray data, Fig. in- 
corporates our whole sample, including the two early-type galax- 
ies. Correspondingly, we have a much better coverage of the low 
mass end of the sample, and so have used a narrower temperature 
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Figure 7. Spatial variation of cumulative total density - normalized to 
the critical density of the Universe - with scaled radius, for the full sam- 
ple, grouped by system temperature. The solid line represents the coolest 
systems (including the two galaxies) (0.3-1.3 keV), increasing in tempera- 
ture through dashed (1.3-2.9 keV), dotted (2.9^4.6 keV), dot-dashed (4.6- 
8 keV) and finally dot-dot-dot-dashed (8-17 keV). 



range for the coolest band (0.3-1.3 keV, instead of 0.3-2.0 keV) to 
provide greater sensitivity to any possible enhancements in density 
concentration at this mass scale. 

The dark matter accounts for the majority of the gravitating 
mass, so it is not surprising that the trend seen in Fig.|6|is present 
in Fig. The difference between the dark matter and total mass 
density seen in the second coolest bin can be attributed to the 
slightly different temperature ranges in each case, used for reasons 
explained previously. 

Although the total density profiles in Fig. Q are essen- 
tially self-similar in the outer regions, there is still a rea- 
sonable degree of scatter. We have compared -R200 with two 
different radii of overdensity, for each object in our X-ray 
sample: we have selected -R2 500 and T^sqq, which are of- 
ten used in the literature (e.g. |A llen 1 Sch midt & Fabianl 1200 it 
iFinoguenov. Reiprich & BohringeJ l200ll) . since X-rav haloes are 
more readily traceable out to these radii than to -R200 ■ We find that, 
on average, an overdensity of 2500 and 500 correspond to 31 and 
66 per cent of -R200, respectively, with corresponding standard de- 
viations across our sample of 5 and 4 per cent. 



5.3 Gas distribution 

Fig. IU shows the gas density profiles for the whole X-ray sample, 
averaged according to their mean X-ray temperature. It is very clear 
that there is a great deal of variation in the five different curves. 
Even for the three profiles representing clusters hotter than 2.9 keV 
there is a variation in the normalization of the lines, although the 
shape of the profiles is very similar. 

The two coolest temperature bands exhibit very different be- 
haviour, having the lowest gas density within 0.3i?200 and show- 
ing no evidence of the central core seen in the hotter clusters. The 
normalization of the cool group lines is consistent with the trend 
towards a lowering of gas density compared to the clusters. This 
behaviour mirrors the trend seen in the gas fraction in Paper 1 and 
provides strong evidence of a deviation from self-similar scaling of 
the IGM both within as well as between haloes of different masses. 



© 2003 RAS, MNRAS OOO.ITWmI 



10 A. J. R. Sanderson and T. J. Ponman 




Temperature (keV) Temperature (keV) 

Figure 9. The ratio of r c to R200 as a function of temperature. The left panel shows all the individual points; the right panel shows the data grouped to a 
minimum of four points per bin. The axes in each plot have been scaled identically. The 1 8 systems with potentially unreliable core radii have been excluded 
(see text for details). 




0.01 0.1 1 

Scaled Radius (R/R 2Q0 ) 



Figure 8. The variation of gas density with scaled radius, for the full sam- 
ple, grouped by system temperature. The solid line represents the coolest 
systems (including the two galaxies) (0.3—1.3 keV), increasing in tempera- 
ture through dashed (1.3-2.9 keV), dotted (2.9^1.6 keV), dot-dashed (4.6- 
8 keV) and finally dot-dot-dot-dashed (8-17 keV). 

5.4 Central density concentration 

Although the underlying gravitational potential in virialized sys- 
tems is expected to be self similar, simulations indicate that the 
concentration of the dark matter should vary slightly with mass 
lNavarroetalll997l) . It was shown in Figs.|6|and0that the coolest 
systems in our sample do indeed appear to be more centrally con- 
centrated - an effect which is most pronounced in the central ~10 
per cent of R200 • 

To understand what is driving this behaviour, we have exam- 
ined the scaling properties of the constant density c ore in the gas 
distribution, r c , as measured from the X-ray data in IPaper it The 
left panel of Fig. [5] shows r c as a fraction of -R200 plotted against 
system temperature. For the hottest clusters the points scatter about 
a self-similar mean value of roughly 10 per cent. As the tempera- 
ture decreases, however, the scatter increases and < 1 keV there is a 
very sharp drop. This trend is in excellent agreement with the pre- 
dictions of the galaxy formation-regulated gas evolution model of 



iBrvanl feOOOl) . as studied in detail bv lWu & Xuel <2002t) : there is a 
remarkable similarity between Fig.|5|and the model and data points 
in fig. 8 from the latter paper. 



To minimize the possibility of systematic bias in this result, we 
have identified and excluded from Fig. [5] those systems for which 
the core radius may be unreliable, owing to the presence of a central 
cooling emission excess. In 17 cases, r c was measured to be smaller 
than the either the size of the radius used to excise the cooling flow 
or the radius within which a cooling flow component was fitted. In 
addition, the cluster Abell 2218 was excluded, since it was neces- 
sary to fix its core radius at the best-fitting value in o rder to st abilize 
the fitting during the calculation of parameter errors JPaper J) , leav- 
ing a total of 48 systems plotted in Fig. [9] To suppress the scatter 
in the relation, we have grouped the points together to a minimum 
of four points per bin, giving a total of 12 bins. The errors on each 
point are determined from the scatter in the X and Y directions. 
This plot is shown in the right panel of Fig.|9| The underlying trend 
in the data is much clearer - it can be seen that the coolest bin is at 
least 10a lower than the flat relation established by the six hottest 
bins. 



It is difficult to explain this behaviour in terms of a fitting 
bias: had these coolest systems exhibited X-ray core radii consis- 
tent with the cluster trend, i.e. > 3 times larger than observed, they 
would very easily have been detected. One possibility is that the 
values of -R200 are anomalously large. To test this we have plotted 
-R200 against temperature in Fig. \H] The slope of this relation is 
somewhat steeper than the self-similar prediction of 0.5, but there 
is clearly no evidence of a trend towards an unusually large -R200 
in the cooler systems, which could account for the discontinuity 
observed in r c /i?200- The most clearly discrepant point at the low- 
mass end is the SO galaxy NGC 1553 (the coolest system in our 
sample), but it appears to have a value of -R200 which is lower than 
its temperature would suggest. However, this object has unusual 
properties, which point to an anomalously high temperature, which 
may have been boosted by energy injection from stellar winds (see 
PaDerj|)- 
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5.5 Star formation efficiency and gas loss from haloes 

To address the issue of star formation efficiency, we have exam- 
ined the star-to-baryon ratio as a measure of the effectiveness with 
which gas has been converted into stellar material. The upper left 
panel of Fig. shows this quantity plotted against system tem- 
perature for the optical sample. It can be seen that there is some 
evidence of a negative correlation, which is significant at the 3.1a 
level. However, this trend is potentially misleading, since it may be 
due at least in part to the variation in gas mass with temperature 
seen above. This means that some allowance must be made for the 
possibility that gas may be lost to the system beyond the virial ra- 
dius, i.e. the gas mass may be reduced without any corresponding 
increase in the stellar mass. 

A better discriminator of star formation efficiency is obtained 
if the stellar mass is normalized to the mass of dark matter, since 
this component is largely immune to bias from non-gravitational 
processes. We have plotted this ratio against X-ray temperature in 
the upper right panel of Fig. 1101 It is clear that the trend seen in 
upper left panel of Fig. llOl has largely vanished, leaving only 0.8cr 
evidence for a systematic variation in star formation efficiency with 
halo temperature. 

To identify the cause of the trend in the top left panel of 
Fig. 1101 we have also plotted the ratio of gas to dark matter mass, 
shown in the lower left panel of Fig.^| A positive correlation be- 
tween this ratio and the X-ray temperature is significant at the 3.3a 
level. This confirms that the variation in gas fraction with tempera- 
ture does indeed fully account for the trend found in the stellar-to- 
baryon ratio. 

5.6 Baryon fraction and constraints on ft m 

Using our measured stellar and gas mass data, we are able to de- 
termine the fraction of mass in baryons, /b, for our optical sample, 
which we show plotted against X-ray temperature in the lower right 
panel of Fig. ^| There is evidence of a modest trend in the data 
(2.5ct significance, using Kendall's K statistic - a non-parametric 
correlation test). The best-fitting power law is represented by the 
solid line, and is given by log / b = (0.37±0.15) log fcT- (0.99± 
0.10); the median value of 0.161 is indicated by the dotted line. 

If virialized systems constitute a fair sample of the baryon con- 
tent of the Universe, we can use this median value to constrain the 
total mass density, Q m = Qb/fb ■ Assuming a baryon density of fib 
= 0.044 ± 0.004 hja~ 2 , from measurements of the po wer spectrum 
of the cosmic microwave background, using WMAP iBennett et alJ 
120031) . we in fer a value of O m = .27 ft 7n . This agrees well with 
the results of lHradeckv et alJfeOOol) . 

If we improve our statistics by combining our mean stellar 
fraction, of (0.034 ± 0.006), with our mean gas fraction for the 
full X-ray sample jPaper J) , of (0.134 ± 0.01) ft 70 3/2 , / b increases 
to 0.168 and f2 m drops to 0.26. Assuming an unbiased measure- 
ment of the stellar and gas fraction, this represents a upper limit 
on f2 m , since any baryonic component to the dark matter would in- 
crease /b. Both these values are in good agreement with the latest 
measu rements of fi m = 0.27 ± 0.04 from WMAP jBennett etail 
l2003t) . 

This li mit on O m i s slightly lower than the value of 0.301q'q3 
inferred by Allen et al. 1 2002), based on a Chandra analysis of six 
massive lensing clusters. The discrepa ncy arises p artly from a dif- 
ference in gas fraction and partly from lAllen et alt s choice of stel- 
lar t o gas ratio. They use a value of 0.19 h ' 5 for the la tter quan- 
tity iFukugita. Hogan & Peebles! 1998UWhite et alll993h . which is 



rather less than our measured value of 0.30 ± 0.04. However, the 
latter value of 0.30 is heavily biased by a handful of groups, which 
have a particularly high fraction of stellar baryons (the uppermost 
points in the top left panel of Fig. fTol : the median stellar to gas ra- 
tio for our sample is 0.21, in much better agreement with existing 
measurements, whic h are based on more massive clusters. 

Moreover, the lAllen et all mean measured gas fraction is 
0.113 ± 0.005, co mpared to our mean of 0.134 ± 0.01 for the 
whole X-ray sample (Paper I). The cause of the discrepancy in both 
cases is the scaling behaviour of the gas fractio n - we find e vidence 
of a rise in gas fraction with radius, whereas lAllen et al] report a 
substantially flat / gas profile in a number of their clusters, albeit 
restricted to a radius of overdensity of ifesoo- By extrapolating a 
constant value out to R200, no allowance is made for an increase 
in /gas, as is predicted by numerical simulations, even when the 
effects of preheating a nd radiative cooling are absent (E ke et alJ 
ll998llFrenk et alll999h . 

Furthermore, we find evidence of a decrease in / gas in cooler 
systems, which has the effect of lowering our mean value. Given the 
possible impact of non-gravitational heating on low mass systems, 
the gas fraction in richer clusters is likely to be a better indicator of 
the universal value. If we calculate an average gas fraction for our 
hottest clusters (>5 keV), we obtain an even higher value, of 0. 17± 
0.01. Combined with our mean stellar fraction from above, this 
places a more stringent upper limit on the mass density, of Q m ^ 
0.22 h^g, which lies just below the WMAP la confidence interval. 



6 DISCUSSION 

6.1 Implications for heating/cooling 

Our results confirm the systematic breaking of self-sim ilarity in the 
IGM that were observed in the gas fraction in Paper I. Clearly the 
thermal history of the hot gas has been altered by the influence 
of non-gravitational physics. The most promising ca ndidates are 
radiative cooling ( Brv arl20 0(tlMuanwong et all200lh and energy 
injection by heating (e.g. Valageas & Silk 1999). Both mechanisms 
are able to account for the X-ray observations: more efficient cool- 
ing in denser (smaller) haloes leads to a depletion of gas in the in- 
ner regions; the energetically boosted IGM is only weakly captured 
in the shallower potential wells of less massive virialized systems, 
thus reducing their gas fraction. 

The apparent self-similarity of the stellar distribution with re- 
spect to the dark matter suggests that star formation has similar 
efficiency in both groups and clusters. However, the behaviour of 
the gas to dark matter ratio points to gas loss in the cooler sys- 
tems, albeit subject to significant extrapolation of the data. To- 
gether, these results are more suggestive of non-gravitational heat- 
ing as the likely mechanism responsible for the observed breaking 
of self-similarity in virialized systems, since this would more natu- 
rally account for depletion of gas in cooler systems without a corre- 
sponding enhanced stellar mass. The alternative, radiative cooling, 
would only be able to reduce gas mass by associated star forma- 
tion, unless it is able to form baryonic dark matter - for example, 
molecular clouds (e.g. Pfenniger, Combes & Martinet 1994). 

The cooling hypothesis (c.f. lKTi^h^^T^omrianin997l) has re- 
ceived increased attention lately, particularly in relation to the ob- 
served entropy proper ties of virialized systems - in deed our own 
results (see Ponman. Sa nderson & Finoguenovl200 3. paper III) are 
consistent with some of the predictions of such models, when the 
effects of feedback from associated star formation are allowed 
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Figure 10. Top left panel: The fraction of baryons in the form of stars, within R200, as a function of system temperature, for the optical sample. Top right: The 
ratio of stellar to dark matter mass within i?200 > as a function of system temperature, for the optical sample. Bottom left: The ratio of gas to dark matter mass 
within R200, as a function of system temperature, for the optical sample. Bottom right: Baryon fraction within R200 as a function of system temperature, for 
the optical sample. The solid line represents the best-fitting power law, which has a logarithmic slope of 0.37 ± 0.15. The dotted line marks the median value 
of 0.161. 



for iVoit & Brvadl200lh . However, the natural consequence of in- 
creased gas cooling is enhanced star formation, unless a substantial 
reservoir of baryons e xists in the f orm of less luminous matter - 
e.g. molecular clouds (Edge 2001). The stellar properties of our 
sample of virialized systems appear consistent with the predictions 
of self-similar scaling, and we find no evidence for a significant 
variation in star formation efficiency across a wide range of halo 
masses. Therefore, we conclude that heating must have had at least 
some role to play in systematically modifying the properties of the 
IGM. 

6.2 Halo formation epoch 

From both the dark matter and total mass density profiles, it is clear 
that there is an enhancement in the central concentration in the very 
coolest systems in our sample. This behav iour is consistent w ith 
hierarchical structure formation (e.g. Blumenth al et alJll984f) . in 
which the smallest haloes collapse at the earliest epochs: the higher 
density of the Universe at this time results in a higher central den- 
sity. Such a systematic variation in the epoch of formation of viri- 
alized systems has previously been observed ISato et all2 000). 



Another possibility is that the greater ages for these sys- 
tems allow more time for accumulat ion of extra ma terial on to 
their haloes, via accretion JSalvador-Sole et alj r"l998X This pro- 
cess would have the effect of increasing the virial radius, but with- 
out significantly modifying the characteristic turnover radius in the 
gas density profile, i.e. the X-ray core radius as a fraction of -R200 
would shrink, as observed in section l5l4l Although small haloes are 
still able to form at the present epoch, there may be an intrinsic bias 
towards observing older, more relaxed systems, which are likely to 
be brighter - indeed, we selected our X-ray sample principally on 
the basis of a regular X-ray morphology. Thus we might expect to 
find fewer examples of cool systems with large a core, compared 
to their virial radius. This may account for the absence of points 
around the cluster trend of r c //?200 ~ 0.1below~l keVinFig.|9| 

The formation epoch of virialized haloes may play quite a 
significant role in influencing their scaling properties. This is par- 
ticular true for nearby systems, where the redshift of observation, 
z bs, is a poor and systematically biased estimator of the redshift 
of formation, Zf . The trend for less massive haloes to be older may 
account for the observed steepening of the M — Tx relation (see 
IPaper J) , in contrast to the apparently self-similar slope found in 
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Figure 11. The variation of -R200 as a function of temperature. The solid 
line depicts the best fitting linear relation, with a slope of 0.60 ± 0.03. 

more distant, massive clusters I Allen etalll200lh . whe re z b B is a 
much better measure of Z{ . However, the simulations of Mathiesen 
1200 lh indicate that this process has a negligible effect on the tem- 
peratures of those clusters which have assembled 75 per cent of 
their final mass by z = 0.6. 



~1 keV. Surprisingly, we find that there is very little variation in 
the shape and normalization of the stellar density profile across the 
sample. However, the gas density clearly departs from this trend: 
the IGM is significantly more extended and less dense in smaller 
haloes. 

We measure a mean stellar mass fraction of (0.032 ± 
0.004) and a median baryon fraction, for our optical sample, of 

— 3/2 

0.161 h 70 . This allows us to place an upper limit of the mass den- 
sity of the Universe of fl m ^ 0.27 hj Q , though it is possible that 
somewhat lower values are indicated by our data: if we restrict our 
analysis to clusters hotter than 5 keV we deduce Q m ^ 0.22 /if . 
This change is a consequence of a trend in / gas with tempera- 
ture, resulting from the influence of non-gravitational processes on 
the intracluster gas. Our results favour energy injection by non- 
gravitational heating as a contributory mechanism for explaining 
the observed breaking of self-similarity in the IGM, since radiative 
cooling alone would lead to an significant increase in star formation 
efficiency in groups, in contrast to our findings. 

As improved X-ray mass estimates are obtained, with XMM- 
Newton and Chandra, and the technology of wide-field optical 
CCD photometry continues to advance, the prospects for very ac- 
curate determination of the mass-to-light ratio and baryon content 
in virialized systems look good. 



7 CONCLUSIONS 

We have conducted a detailed study of the mass composition of a 
large sample of virialized systems. Using X-ray data for our whole 
sample of 66 objects, we have derived the gas and total gravitating 
mass profiles as a function of radius. Optical data for a subsample 
of 32 groups and clusters have allowed us to determine the stellar 
mass distribution, and thus infer the dark matter density profile. 

We have determined the deprojected luminosity distribution 
from unweighted surface density fits to galaxy positions from the 
APM survey, combined with aperture luminosity measurements 
taken from a small number of sources in the literature. Using the 
galaxy density profile we are able to extrapolate the light to our 
nominal virial radius of -R200 (as determined from the X-ray data) 
to yield the deprojected optical luminosity, in the Bj photometric 
band. 

The scaling properties of the total B band optical luminosity 
within -R200 are consistent with self-similarity; M oc [J- 08 ± - 12 
and L oc y 1(32 ± 014 ; f or t ( a ] m ass and emission-weighted tem- 
perature, respectively. Similarly, the mass-to-light ratio remains es- 
sentially constant across the sample, showing no trend with X-ray 
temperature (M/L oc r p-°- 0fi±0 - 17 y^ we g n( j a logarithmic mean 
M/L of 233^28 ^ 7 o (^/^B.j)©, in good agreement with other 
measurements based on X-ray mass estimates. Our data therefore 
provide no evidence for a significant increase in star formation rate 
in galaxy groups, which is constant across most of the sample. 

We find that the dark matter is the most centrally concentrated 
mass component, followed by the galaxies, with the X-ray emitting 
IGM the most extended component. Comparing the spatial distri- 
bution of these mass components across two decades of halo mass, 
we find that the dark matter and total mass density profiles of our 
sample are nearly self-similar, but for a clear central excess in the 
coolest systems (< 1.5 keV). We attribute this enhancement in cen- 
tral density concentration to a sharp decline in the size of the gas 
core radius, r c (normalized to -R200) amongst systems cooler than 
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